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Abstract 

We investigate the QCD chiral phase transition at high baryon number den- 
sity within the linear quark meson model for two flavors. The method we em- 
ploy is based on an exact renormalization group equation for the free energy. 
Truncated nonperturbative flow equations are derived at nonzero chemical 
potential and temperature. Whereas the renormalization group flow leads to 
spontaneous chiral symmetry breaking in vacuum, we find a chiral symme- 
try restoring first order transition at high density. Combined with previous 
investigations at nonzero temperature, the result implies the presence of a 
tricritical point with long-range correlations in the phase diagram. 
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I. INTRODUCTION 



The behavior of QCD at high temperature and baryon density is of fundamental interest 
and has applications in cosmology, the astrophysics of neutron stars and the physics of 
heavy ion collisions. Over the past years, considerable progress has been achieved in our 
understanding of high temperature QCD, where simulations on the lattice and universality 
arguments played an essential role. Recently, nonperturbative flow equations, based on the 
Wilsonian formulation of the renormalization group, were derived for an effective linear quark 
meson model |]-f3|. The results of this approach account for both, the low temperature chiral 
perturbation theory domain of validity as well as the high temperature domain of critical 
phenomena for two quark flavors 0. The method may help to shed some light on the 
remaining pressing questions at high temperature, like the nature of the phase transition for 
realistic values of the strange quark mass. 

Our knowledge of the high density properties of strongly interacting matter is rudimen- 
tary so far. There are severe problems to use standard simulation algorithms at nonzero 
chemical potential on the lattice because of a complex fermion determinant. Different non- 
perturbative methods as the Wilsonian "exact renormalization group" or Schwinger-Dyson 
equations seem to present promising alternatives. It is the purpose of this letter to pursue 
the former approach based on an exact nonperturbative flow equation for a scale dependent 
effective action which is the generating functional of the 1PI Green functions in 

the presence of an infrared momentum cutoff ~ k. The renormalization group flow for the 
average action interpolates between a given short distance or classical action S and the 
standard effective action T, which is obtained by following the flow for r& to k = 0, thus 
removing the infrared cutoff in the end. 

QCD at nonzero baryon density is expected to have a rich phase structure with different 
possible phase transitions as the density is increased. A prominent example is the liquid-gas 
nuclear transition which has been studied in low energy heavy ion collisions. Large efforts in 
ultra-relativistic heavy ion collision experiments focus on a transition at high temperatures 
and / or densities into what is generally known as the quark-gluon plasma - a new phase of 
matter in which color is screened rather than confined and chiral symmetry is restored. In 
addition to the nuclear and quark matter phases a number of possibilities like the formation 
of meson condensates or strange quark matter have been investigated. Recently, an extensive 
discussion has focused on the symmetry of the high density ground state, where condensates 
of quark Cooper pairs may break the color symmetry of QCD spontaneously and lead to a 
superconducting phase [|],[5| . We will concentrate in this work on the chiral phase transition 
at nonzero baryon number density and apply our methods to an effective linear quark meson 
model for the low energy degrees of freedom of QCDQ We note that the methods we discuss 
are equally suited to treat quark pair condensation and leave this for future work|. 



First calculations indicate that diquark condensates have only little influence on the order pa- 
rameter associated with chiral symmetry breaking in the vacuum ||. 

2 In the context of quark pair formation the renormalization of quark operators in the vicinity of 
the Fermi surface has been studied in ||. 
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The linear quark meson model we consider can be viewed as a generalization of the 
Nambu-Jona-Lasinio model for QCD 0. Though these approaches do not explain the con- 
finement of quarks, they are adequate for many purposes and their predictions show convinc- 
ing agreement with the results of real and numerical experiments. Calculations at nonzero 
baryon density are typically based on a mean field approximation and it was claimed |7[] that 
the order of the chiral phase transition is ambiguous. Recently, the order of this high density 
transition has raised considerable interest. One can argue on general grounds that if 
two-flavor QCD has a second order transition (crossover) at high temperatures and a first 
order transition at high densities, then there exists a tricritical point (critical endpoint) with 
long-range correlations in the phase diagram. The physics around this point is governed by 
universality and may allow for distinctive signatures in heavy ion collisions [f||8],|9|] . 

In this work it is our main concern to show that in a proper treatment beyond mean 
field theory the order of the chiral phase transition can be fixed within the models under 
investigation. For this purpose, a crucial observation is the strong attraction of the flow 
to partial infrared fixed points @||. The two remaining relevant or marginal parameters 
can be fixed by the phenomenological values of f„ and the constituent quark mass. For 
two massless quark flavors we find that chiral symmetry restoration occurs via a first order 
transition between a phase with low baryon density and a high density phase. 

We emphasize, nevertheless, that the linear quark meson model captures the low density 
properties of QCD only incompletely since the effects of confinement are not included. In 
particular, for a discussion of the liquid-gas nuclear transition the inclusion of nucleon 
degrees of freedom seems mandatory ||10|| . On the other hand this model is expected to 
provide a reasonable description of the high density properties of QCD. 



II. LINEAR QUARK MESON MODEL 

It is an important property of the average action r& that at any scale k only fluctuations 
with momenta in a small range around k influence its flow with k. This implies that for each 
range of k only those degrees of freedom have to be included which are relevant around k. In 
QCD these may comprise compounds of quarks and gluons. As long as one concentrates only 
on the chiral properties of QCD one may employ for scales k below a "compositeness scale" 

a description in terms of quark and (pseudo-) scalar meson degrees of freedom [|TT1,p|,|3l| . 
All other degrees of freedom are assumed to be integrated out. The scalar and pseudoscalar 
mesons are described by a complex field $ which transforms as (N, N) under the chiral 
flavor group SUl{N) x SUr(N) with N the number of light quark flavors. Chiral symmetry 
breaking to a vector-like subgroup SUv{N) occurs through a nonvanishing chiral condensate 
= a S ab . 

We will restrict our discussion here to two massless quarks. A more realistic treatment 
would have to include the small up and down quark masses and the finite, but much heavier, 
strange quark mass. Though their inclusion will change certain quantitative estimates, they 
most likely do not change the qualitative outcome of the investigation discussed below. In 
the vacuum the scalar iso-triplet as well as the pseudoscalar singlet (associated with the rj' 
meson) contained in <3> are significantly heavier than the pions and the sigma meson. This 
property even holds in the symmetric phase where the pions are no longer Goldstone bosons. 
We take this as a guide also at nonzero density and decouple these degrees of freedom. For 
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N = 2 this can be achieved in a chirally invariant way. This results in a parameterization 
of $ in terms of the cr-meson field and the three pion field components tx 1 as 



(1) 



where n (Z = 1, 2, 3) denote the Pauli matrices. 

Closed nonperturbative flow equations follow from the exact renormalization group equa- 
tion by a suitable truncation of the effective average action. We consider here a rather 
simple truncation of the momentum dependence in terms of standard kinetic terms pa- 
rameterized by a running meson wave function renormalization constant Z^^ and a scale 
dependent Yukawa coupling h^. For the meson interaction we consider the most general 
form of the scalar potential Uj. consistent with the symmetries. For two flavors chiral sym- 
metry implies that Uk depends on only one invariant p = Tr($^$). The effective potential 
U = Um.k->oUk = linifc^o Tfc T/V is the relevant quantity for a study of the spontaneous 
breaking of chiral symmetry and is obtained for constant field configurations. At its mini- 
mum p the effective potential is related to the energy density e, the entropy density s, the 
quark number density n q and the pressure p by 



U(p ;fx,T) 



Ts 



jjn q 



-P 



(2) 



Our truncation precludes any nontrivial momentum dependence of the effective quark and 
meson propagators or the interactions. In particular, it clearly does not account for con- 
finement effects. At nonzero temperature T and chemical potential /1 associated to the 
conserved quark number our ansatz for r*. reads 



l/T 



dx° / d 3 x\ ^V<^ + P7>a + h k ij> 



1 + 7 5 ». 1 - 7 5 + 



(3) 



We note that in the Euclidean formalism nonzero temperature results in (anti-)periodic 
boundary conditions for (fermionic) bosonic fields in the Euclidean time direction with peri- 
odicity l/T. A nonzero chemical potential \i to lowest order results in the term ~ ipijj ~f°ijj a 
appearing on the right hand side of (|3|). Our approximation neglects the dependence of Z$ fc 
and hk on p and T. We also neglect a possible difference in normalization of the quark kinetic 
term and the baryon number current. The form of the effective action at the compositeness 
scale, $], serves as an initial value for the renormalization group flow of r^^, $] for 

k < We will consider here the case that Z^^ "C 1. The limiting case Zq,^ = can be 
considered as a solution of the corresponding Nambu-Jona-Lasinio model where the effective 
four-fermion interaction has been eliminated by the introduction of auxiliary meson fields. 



III. FLOW EQUATION FOR THE EFFECTIVE POTENTIAL 



The dependence on the infrared cutoff scale k of the effective action r*. is given by an 
exact flow equation |I,12|, which for fermionic fields ip (quarks) and bosonic fields $ (mesons) 
reads 
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(4) 



(2) 

Here 1^ is the matrix of second functional derivatives of with respect to both fermionic 
and bosonic field components. The first trace on the right hand side of (^) effectively 
runs only over the bosonic degrees of freedom. It implies a momentum integration and a 
summation over flavor indices. The second trace runs over the fermionic degrees of freedom 
and contains in addition a summation over Dirac and color indices. The exact flow equation 
closely resembles a one-loop equation with the difference that the full inverse propagator 

(2) 

T k at the scale k appears instead of the classical propagator. The infrared cutoff function 
R k has a block substructure with entries RkB and R k p. Here R k B denotes the infrared cutoff 
function for the bosonic fields and we employ an exponential cutoff function 

_ Z^ h q 2 

KkB ~ exp(q 2 /k 2 ) - 1 • (5j 

With this choice, R k s acts as an additional mass term R k g — Z^ k k 2 for the low momentum 
q 2 <C k 2 modes. The infrared cutoff function for the fermions R k p should be consistent with 
chiral symmetries. This can be achieved if R kF has the same Lorentz structure as the kinetic 



term for free fermions In presence of a chemical potential \i we use 

RkF = (Yq^ + i^l Q )r k F ■ (6) 

The effective squared inverse fermionic propagator is then of the form 

PkF = [(go + + <f](l + r kF ) 2 

= (g + i^f + q 2 + k 2 Q{k% - (q + ip) 2 - q 2 ) , (7) 

where the second line defines r k p and one observes that the fermionic infrared cutoff acts as 
an additional mass-like terrn^ ~ k 2 Here the B-function can be thought of as the limit of 
some suitably regularized function, e.g. 6 e = [exp{(g + ip) 2 + q 2 — k%} fe + l] -1 . 

We compute the flow equation for the effective potential U k from equation @ using the 
ansatz (|3|) for Y k and we introduce a renormalized field p = Z$ tk p and Yukawa coupling 

— 1/2 — 

h k = k ' h k . The flow equation for U k 

d d 

— U k (p; T, p) = -^rU kB {p] T, p) + T^U kF (p; T, p) (8) 
obtains contributions from bosonic and fermionic fluctuations, respectively, 



3 The exponential form (|5|) of the cutoff function R k B renders the first term on the right hand side 
of (||) both infrared and ultraviolet finite. No need for an additional ultraviolet regularization arises 
in this case. This is replaced by the necessary specification of an initial value at the scale k&. 
Here /c<j> is associated with a physical ultraviolet cutoff in the sense that effectively all fluctuations 
with q 2 > k\ are already included in T k(j? . A similar property for the fermionic contribution is 
achieved by the 0-function in (fj). 
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d_ 

dk 



U kB {p] T, p) 



I T V f 1 ^W) f 3 

2 (2tt) 3 ^ dk \P kB (q 2 ) + U' k (p;T,p) 



+ 



1 



-^U kF (p;T, p) 



P kB (g 2 ) + U' k ( P] T, p) + 2pU' k \p- T,p) j 

kQ(kl- (q + tp) 2 - q 2 ) 



-8N C TJ2 



(2tt) 3 P fcF ((g + ip) 2 + q 2 ) + Hp/2 



(9) 
(10) 



Here q 2 = q 2 + g 2 with qo(n) = 2mrT for bosons, qo(n) = (2n + 1)ttT for fermions (n G ZZ) 
and A^ c = 3 denotes the number of colors. The scale dependent propagators on the right 
hand side contain the momentum dependent pieces P kB = q 2 + Z^ k R k (q 2 ) and P k F given 
by (^|) as well as mass terms, where U' k , U k denote the first and second derivative of the 
potential with respect to p. The only explicit dependence on the chemical potential p 
appears in the fermionic contribution (|I0| ) to the flow equation for U k . It is instructive to 
perform the summation of the Matsubara modes explicitly for the fermionic part. Since the 
flow equations only involve one momentum integration, standard techniques for one loop 
expressions apply and we find 



d_ 

dk 



U kF (p;T,p) = -8N C 



d 4 q 



k 6(4 



Q 



(2tt) 4 g 2 + k 2 + h 2 k p/2 



4:N r 



d 3 q 



( 27 ^) 3 Jq 2 + k 2 + h 2 k p/2 



x 



cxp 



q 2 + k 2 + h 2 p/2 - p)/T 



1 exp 



q 2 + k 2 + h 2 p/2 + p)/T 



(11) 



where, for simplicity, we sent — >• oo in the p, T dependent second integral. This is justified 
by the fact that in the p, T dependent part the high momentum modes are exponentially 
suppressed. 

For comparison, we note that within the present approach one obtains standard mean 
field theory results for the free energy if the meson fluctuations are neglected, dU kB /dk = 0, 
and the Yukawa coupling is kept constant, h k = h in (|TT1) . The remaining flow equation 
for the fermionic contribution could then easily be integrated with the (mean field) initial 
condition U k<6 (p) = m| p. In the following we will concentrate on the case of vanishing tem- 
perature^]. We find (see below) that a mean field treatment yields relatively good estimates 
only for the /i-dependent part of the free energy U(p; p) — U (p; 0). On the other hand, mean 
field theory does not give a very reliable description of the vacuum properties which are 
important for a determination of the order of the phase transition at p ^ 0. 



For a solution of the flow equations at nonzero temperature and vanishing chemical potential 
sec 1. 
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IV. RENORMALIZATION GROUP FLOW AT NONZERO CHEMICAL 

POTENTIAL 



In the limit of vanishing temperature one expects and observes a non-analytic behavior 
of the /z-dependent integrand of the fermionic contribution ([□]) to the flow equation for U k 
because of the formation of Fermi surfaces. Indeed, the explicit /i-dependence of the flow 
equation reduces to a step function 



d_ 

dk 



■U kF (r,n)=-8N e ] (2^2 + ^2/2 

— OO 

OO _« _ 



+ 4N < I (fli ! k e(fJi-Jq* + k* + hlp/2). (12) 

-L (27T) ^ 2 + k 2 + Hp/2 V v ; 

The quark chemical potential p enters the bosonic part (|9p of the flow equation only im- 
plicitly through the meson mass terms U k (p; p) and U' k (p; p) + 2pUjt(p; p) for the pions and 
the cr-meson, respectively. For scales k > p the ©-function in ([12]) vanishes identically and 
there is no distinction between the vacuum evolution and the p ^ evolution. This is due 
to the fact that our infrared cutoff adds to the effective quark mass (k 2 + h 2 k p/2) 1 / 2 . For a 
chemical potential smaller than this effective mass the "density" —dUk/dp vanishes whereas 
for larger p one can view p = {cfp(p, k, p) + k 2 + ti 2 k p/2] 1 / 2 as an effective Fermi energy for 
given k and p. A small infrared cutoff k removes the fluctuations with momenta in a shell 
close to the physical Fermi surfaceQ/i 2 — h 2 k p/2 — k 2 <q 2 <p 2 — h 2 k=0 p/2. Our flow equation 
realizes the general idea |L3| that for p ^ the lowering of the infrared cutoff k —>■ should 
correspond to an approach to the physical Fermi surface. For a computation of the meson 
effective potential the approach to the Fermi surface in fll2"|) proceeds from below and for 
large k the effects of the Fermi surface are absent. By lowering k one "fills the Fermi sea". 

In vacuum the evolution of Uk, the Yukawa coupling hk and the meson wave function 
renormalization Zq> )k have been computed earlier [@,[3|](]. The flow equations for hk and Zq,^ 
can be derived from (|j) using the truncation (^) and one finds that for large enough k their 
running is well approximated by 0|| {k < k§) 

One observes a very strong Yukawa coupling hk for k in the vicinity of kq>. (The Landau 
pole at k = corresponds to taking —►().) This strong Yukawa coupling implies 



If one neglects the mesonic fluctuations one can perform the fe-integration of the flow equa- 
tion (|l2| ) in the limit of a /^-independent Yukawa coupling. One recovers (for fc| S> k 2 + h 2 p/2, p 2 ) 
mean field theory results except for a shift in the mass, h 2 p/2 — > h 2 p/2 + k 2 , and the fact that 
modes within a shell of three-momenta p 2 — h 2 p/2 — k 2 < q 2 < p 2 — h 2 p/2 are not yet included. 
Because of the mass shift the cutoff k also suppresses the modes with q 2 < k 2 . 

6 In (2],|j a different variant of the fermionic infrared cutoff function RkF than (|7|) was used. 
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a fast approach of running couplings to partial infrared fixed points [@,|J. More precisely, 
for a rescaled potential Uk(p)/k A = J2™=o h u k\p/^ 2 ) n ^ ne running dimensionless couplings 

(n) 

u k for n > 2 approach values proportional to powers of the Yukawa coupling h k , with 
v$ /h\ ~ 1, uf'/hl ~ —0.0164, u k /h\ ~ 0.0105 and similarly for the higher couplings |§. 
In consequence, the details of the short distance meson interactions become unimportant 
and besides the marginal coupling hj. the only relevant parameter of the model corresponds 
to the mass parameter . We note that the large value of h k is phenomenologically 
suggested by the comparably large value of the constituent quark mass. The observed fixed 
point behavior allows us to fix the model by only two phenomenological input parameters 
and we use f n = 92.4 MeV and 300 MeV^M g ^350 MeV. Previous results for the evolution 
in vacuum HH show that for scales not much smaller than fc$ ~ 650 MeV chiral symmetry 
remains unbroken. This holds down to a scale k x sB — 450 MeV at which the meson potential 
Uk{p) develops a minimum at po,fc > thus breaking chiral symmetry spontaneously. Below 
the chiral symmetry breaking scale running couplings are no longer governed by the partial 
fixed point. In particular, for k^k x ss the Yukawa coupling hf. and the meson wave function 
renormalization depend only weakly on k and approach their infrared values. 

At p ^ we will follow the evolution from k = k x $B to k = and neglect the /in- 
dependence of hk and Z<$>^ i n this range. According to the above discussion the initial 
value Uk xSB is //-independent for p < k x sB- At this scale SB is well approximated by 

the polynomial form where the coefficients u ksB , n = 2,3,4, are estimated by their partial 
infrared fixed point values and coefficients with n > 4 are neglected. The scalar mass 
parameter u k crosses from positive to negative values at the chiral symmetry breaking 
scale and thus vanishes at k = k x ss- The constant u k °^ sg is fixed such that the pressure of 
the physical vacuum is zero. We do not assume a polynomial form for U k (p; p) for k < k x sB 
and solve the flow equation (§) with @, fljjD numerically^ as a partial differential equation 
for the potential depending on the two variables p and k for given p. 

In the fermionic part ([12]) of the flow equation the vacuum and the /i-dependent term 
contribute with opposite signs. This cancelation of quark fluctuations with momenta below 
the Fermi surface is crucial for the restoration of chiral symmetry at high density. In vacuum, 
spontaneous chiral symmetry breaking is induced in our model by quark fluctuations which 
drive the scalar mass term U' k (p = 0) from positive to negative values at the scale k = k xSB . 
(Meson fluctuations have the tendency to restore chiral symmetry because of the opposite 



Nonzero current quark masses, which are neglected in our approximation, result in a pion mass 
threshold and would effectively stop the renormalization group flow of renormalized couplings at 
a scale around f%. To mimic this effect one may stop the evolution by hand at kf ~ m n . We 
observe that our results are very insensitive to such a procedure which can be understood from the 
fact that a small infrared cutoff - induced by kf or by the nonzero current quark masses - plays 
only a minor role for a sufficiently strong first order transition. For the results presented in the 
next section we use kf = 100 MeV. The flow of the potential U k around its minima and for its 
outer convex part stabilizes already for k somewhat larger than kf. Fluctuations on larger length 
scales lead to a flattening of the barrier between the minima. The approach to convexity is not 
relevant for the present discussion. 
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FIG. 1. The zero temperature effective potential U (in MeV 4 ) as a function of <r = (p/2) 1 / 2 for 
different chemical potentials. One observes two degenerate minima for a critical chemical potential 
jj, c /M q = 1.025 corresponding to a first order phase transition at which two phases have equal 
pressure and can coexist (M q = 316.2 MeV). 

relative sign, cf. (^).) As the chemical potential becomes larger than the effective mass 
(k 2 + h\p/2) l l 2 quark fluctuations with momenta smaller than q^(/i, k, p) = jj 2 — k 2 — hip/ 2 
are suppressed. Since cfp is monotonically decreasing with p for given p and k the origin of 
the effective potential is particularly affected. We will see in the next section that for large 
enough p this leads to a second minimum of Uk=o(p; p) at p = and a chiral symmetry 
restoring first order transition. 



V. HIGH DENSITY CHIRAL PHASE TRANSITION 



In vacuum or at zero density the effective potential U(a), a = dp/2, has its minimum 
at a nonvanishing value a = f n /2 corresponding to spontaneously broken chiral symmetry. 
As the quark chemical potential p increases, U can develop different local minima. The 
lowest minimum corresponds to the state of lowest free energy and is favored. In Figure [I] 
we plot the free energy as a function of a for different values of the chemical potential 
p = 322.6, 324.0, 325.2 MeV. Here the effective constituent quark mass is M q = 316.2 MeV. 
We observe that for p < M q the potential at its minimum does not change with p. Since 
(cf. (i) 

«, = ~T (") 

|min 

we conclude that the corresponding phase has zero density. In contrast, for a chemical 
potential larger than M q we find a low density phase where chiral symmetry is still broken. 



S 




H/M q 



FIG. 2. The plot shows n q , where n q denotes the quark number density as a function of [i in 
units of the effective constituent quark mass (M q = 316.2 MeV). 

The quark number density as a function of p is shown in Figure |2|. One clearly observes the 
non-analytic behavior at p — M q which denotes the "onset" value for nonzero density. From 
Figure |l| one also notices the appearance of an additional local minimum at the origin of U. 
As the pressure p = —U increases in the low density phase with increasing p, a critical value 
p c is reached at which there are two degenerate potential minima. Before p can increase 
any further the system undergoes a first order phase transition at which two phases have 
equal pressure and can coexist. In the high density phase chiral symmetry is restored as can 
be seen from the vanishing order parameter for p > p c . The curvature of the potential at 
the origin for p^p c is d 2 U/da 2 (a = 0) ~ (1.45M g ) 2 . We note that the relevant scale for 
the first order transition is M q . For this reason we have scaled our results for dimensionful 
quantities in units of M q . 

For the class of quark meson models considered here (with M q j 'f n in a realistic range 
around 3-4) the first order nature of the high density transition has been clearly established. 
In particular, these models comprise the corresponding Nambu-Jona-Lasinio models where 
the effective fermion interaction has been eliminated by the introduction of auxiliary bosonic 
fields. Our treatment, which is based on a nonperturbative evolution equation for the free 
energy, goes beyond mean field theory and puts it into a more systematic context. For the 
quark meson model this method has been used in the past to study the chiral phase transition 
at nonzero temperature, including the critical properties near the second order transition 
for two massless flavors characteristic of the three dimensional 0(4) universality class |§. In 
the chiral limit, we therefore find a second order chiral transition at zero chemical potential 
and a first order chiral transition at zero temperature. By continuity these transitions meet 
at a tricritical point in the (/i,T)-plane Away from the chiral limit, the second order 

chiral transition turns into a smooth crossover. The first order line of transitions at low 
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temperatures now terminates in a critical endpoint with long-range correlations |5]||. Using 
the method employed here a precision estimate of the universal critical equation of state of 



this Ising endpoint has been given in |14 



The extent to which the transition from nuclear matter to quark matter in QCD differs 
from the transition of a quark gas to quark matter in the quark meson model or NJL- 
type models has to be clarified by further investigations. At the phase transition, the quark 
number density in the symmetric phase n 1 ^ = 0.52M g turns out in this model to be not much 
larger than nuclear matter density, n^ uc = 152 MeV. Here the open problems are related to 
the description of the low density phase rather than the high density phase. Quite generally, 
the inclusion of nucleon degrees of freedom for the description of the low density phase shifts 
the transition to a larger chemical potential and larger baryon number density for both the 
nuclear and the quark matter phases |)10| |. An increase in fi c results also from the inclusion of 



vector mesons. This would not change the topology of the phase diagram inferred from our 
model. Nevertheless, the topology of the QCD phase diagram is also closely connected with 
the still unsettled question of the order of the high temperature (/i = 0) transition for three 
flavors with realistic quark masses. If the strange quark mass is too small, or if the axial 
U(l) symmetry is effectively restored in the vicinity of the transition, then one may have 
a first order transition at high T which is driven by fluctuations [HJ. In this case, critical 



endpoints only occur if the two first order regions atT>// and /i ^> T are disconnected. 

Finally, the low-density first order transition from a gas of nucleons or the vacuum to 
nuclear matter (nuclear gas-liquid transition) can only be understood if the low-momentum 
fermionic degrees of freedom are described by nucleons rather that quarks |ID|. The low- 



density branch of Figure |2| cannot be carried over to QCD. We emphasize that nucleon 
degrees of freedom can be included in the framework of nonperturbative flow equations. 
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